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Density functional theory (DFT) is an incredible success story. The low computational cost, 
combined with useful (but not yet chemical) accuracy, has made DFT a standard technique in most 
branches of chemistry and materials science. Electronic structure problems in a dazzling variety of 
fields are currently being tackled. However, DFT has many limitations in its present form: Too many 
approximations, failures for strongly correlated systems, too slow for liquids, etc. This perspective 
reviews some recent progress and ongoing challenges. 



I. INTRODUCTION 

Over the past 20 years, DFT has become a much used 
tool in most branches of chemistry. Many experimen- 
tal investigations in organic and inorganic chemistry rou- 
tinely include such calculations, using a popular code, a 
standard basis, and a standard functional approximation 
PP. A similar transformation is now underway in ma- 
terials science where, in the past decade, improvements 
in both hardware and codes have made it possible to 
perform systematic comparisons with experiment across 
large ranges of materials, learning which approximations 
work and why, and allowing for true first-principles pre- 
dictions of properties. Among notable recent successes 
are the prediction of new catalysts [2] and new Li bat- 
tery materials [3] in the Materials Genome Project. A 
complementary aspect of this story is shown in Fig. [T] 
which plots the number of papers given by Web of Knowl- 
edge when DFT is searched as a topic (grey bars). This 
will soon reach 10,000 per year, vindicating the 1998 No- 
bel prize in chemisty, which went to Walter Kohn [3] for 
inventing the theory and to John Pople [5] for making it 
accessible through popular computational packages. The 
figure also marks the fraction of papers citing B3LYP 6 - 
|S], currently the most popular approximation in chem- 
istry, and PBE [9], the most popular approximation in 
materials. Clearly, applications to materials will soon 
outstrip those in chemistry. 

This perspective is for a general audience, and focuses 
on fundamental general aspects of DFT, rather than de- 
tailed computational procedures and results for specific 
systems. Because DFT is now applied so broadly, no 
such article can hope to be comprehensive. The topics 
covered here are designed to give a flavor of how the field 
works, and are mostly those I have personally worked 
in. Great longer reviews exist for experts [I0HT2] and for 
users [T5]. as well as introductions for neophytes [TH ITS] . 
An excellent complementary perspective, including com- 
putational issues and non-DFT approaches, is Ref. [Ttj] . 
I work always within the Born-Oppenheimer approxima- 
tion, the non-relativistic limit, and discuss only the basis- 
set limit of calculations. I use atomic units (lengths in 
Bohr, energies in Hartree) except where otherwise speci- 
fied. For simplicity, I refer always to density functionals, 
but modern calculations use spin-density functional the- 
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FIG. 1. Numbers of papers when DFT is searched as a topic 
in Web of Knowledge (grey), B3LYP citations (blue), and 
PBE citations (green, on top of blue). 



ory, a generalization [17] . 



II. A BRIEF HISTORY 

Our story begins in 1926 with the creation of Thomas- 
Fermi theory [181 119) , an approximate method for find- 
ing the electronic structure of atoms using just the one- 
electron ground-state density, p(r), but too crude to bind 
molecules [20] . In the 50's, Slater [21] intuitively com- 
bined this idea with Hartree's orbital method [25] in the 
Xcv scheme. Later, the Hohenberg-Kohn (HK) theorem 
[2"5] proved that an exact method based on p(r) exists 
in principle. The modern version in use today is Kohn- 
Sham (KS) DFT, which defines self-consistent equations 
that must be solved for a set of orbitals whose density, 
p(r) is defined to be exactly that of the real system [2"4"]. 
In these equations, a small but vital contribution to the 
energy, the exchange-correlation (XC) energy, must be 
given in terms of p(r). In principle, and for small systems, 
this functional can be found exactly, but turns out to be 
more expensive than direct solution of the Schrbdinger 
equation [25] . In practical calculations, the XC contribu- 
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tion is approximated, and the results are only as good as 
the approximation used. 

The simplest XC approximation is the local density ap- 
proximation (LDA) 2 1 which became the popular stan- 
dard in calculations on solids in the 70's and 80's. But 
molecules in LDA are typically overbound by about 1 
eV/bond, and in the late 1980's generalized gradient ap- 
proximations (GGA's) [IB] produced an accuracy that 
was useful in chemical calculations. In the early 90's, hy- 
brids were introduced by Becke [B], replacing a fraction 
of GGA exchangejT] with Hartree-Fock exchange, leading 
to the ubiquitous B3LYP [SJ, the most popular approx- 
imation in use in chemistry today. On the other hand, 
the PBE GGA [3] has come to dominate applications 
to extended systems (materials). We denote these three 
(LDA, PBE, and B3LYP) as the standard approxima- 
tions, meaning that they are the most popular examples 
of each type of approximation, and dominate the user 
market, as shown by Fig. [l] 

To give an idea of how much (or how little) progress 
is made in DFT development, we go back to the year 
of the Nobel prize, 1998. Some of the most prominent 
practical difficulties from back then include (i) gaps of 
bulk solids are underestimated, (ii) van der Waals miss- 
ing from popular functionals, (iii) strongly correlated sys- 
tems poorly treated, and (iv) no good scheme for excita- 
tions. Of course, there are many others, but this sample 
will give us some idea of how things work. The rest of 
this article interweaves sections labeled progress, indicat- 
ing areas where substantial progress has been made since 
that time, with general background, explaining relevant 
concepts, and challenges, areas that need improvement 
and where we can hope for progress in the next decade. 



III. PROGRESS: MATERIALS AND 
NANOSCIENCE 

As shown in Fig. [T] materials applications now share 
the limelight with chemistry. Nanoscience completely 
interweaves these two, such as when a molecule is ad- 
sorbed on a surface. Many areas of materials research, 
especially those related to energy, desperately need input 
from electronic structure methods. Most calculations of 
materials use codes that differ from traditional quantum 
chemistry codes, because they employ plane waves that 
satisfy periodic boundary conditions [37] . Such calcu- 
lations converge much faster to the bulk limit than by 
taking ever larger finite clusters of the material. In fact, 
almost all popular codes are designed either for finite 
molecular systems [25J or for extended bulk systems [2H] , 
although codes that treat both are beginning to appear 

PUGS]. 

In solid-state physics, DFT has always been more pop- 
ular than traditional direct solution of the Schrodinger 
equation, because the Hartree-Fock approximation has 
unpleasant singularities for zero-gap materials, i.e., met- 
als. The next logical step [32] beyond the standard 



approximations are meta-GGA's, which include the ki- 
netic energy density as an input, and can yield accu- 
rate ground-state energies simultaneously for molecules, 
solids, and surfaces [33J. But an outstanding failure has 
been DFT's inability to provide good estimates of the 
fundamental gaps of semiconductors and insulators, a 
crucial quantity for much materials research, such as im- 
purity levels in doped semiconductors |34j . 

The fundamental (or transport) gap is / — A, where 
/ is the ionization potential and A the electron affinity. 
It is well-established that the KS gap, i.e., the differ- 
ence between the KS HOMO and LUMO energies, does 
not equal this value, even with the exact XC functional 
[551 13B] . Calculations suggest it is typically substantially 
smaller (by about 50%) [37J. The LDA and GGA ap- 
proximations yield accurate values for the KS gap. But 
if / — A were calculated by adding and subtracting elec- 
trons from a large cluster of the material, it would match 
the KS gap within such approximations, because they 
incorrectly allow electrons to completely delocalize over 
insulating solids. Such approximations lack a derivative 
discontinuity that accounts for this difference |36| . 

A practical way around this problem is to use a hy- 
brid, but treating the exchange term as orbital depen- 
dent, via generalized KS theory [35] . With mixing pa- 
rameters of about 20%, this typically yields much better 
gaps for semiconductor materials. A recent functional, 
called HSE [3j|l HO], not only mixes in some HF, but also 
performs a range-separation. Based on exact theorems 
of SavinJJT], the short-range part of the HF is treated 
exactly, while the long-range contribution is treated by 
approximate DFT. The resulting functional, used to cal- 
culate gaps in the generalized scheme, appears to work 
accurately for a large variety of moderate-gap semicon- 
ductors [40], overcoming the problems of LDA and GGA. 
It yields accurate fundamental gaps when excitonic ef- 
fects are negligible, and is closer to optical gaps when 
they are not [421 143]. 

This kind of progress is very welcome in an era in 
which we must accurately tackle moderately correlated 
systems. Transition metal oxides play vital roles in many 
energy- related problems, such as creation of efficient 
photovoltaics [31] . The standard approximations over- 
delocalize the d-electrons, leading to highly incorrect de- 
scriptions. Many practical schemes (HSE 39 , GGA+U 
03 US], Dynamical Mean-Field Theory (DMFT) [47] ) 
can correct these difficulties, but none has yet become 
a universal tool of known performance for such sys- 
tems. Very recently, a promising non-empirical scheme 
has been suggested for extracting gaps using any approx- 
imate functional [35] . 

An alternative approach to directly tackling such prob- 
lems is to study them in simpler situations, and test 
suggested remedies on cases where exact, unambigu- 
ous answers can be easily obtained. Recently, an ex- 
tremely powerful technique for direct solution of many- 
body problems, called DMRG, has been adapted to tackle 
a one-dimensional analog of the real world [4jJ] . DMRG 
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is powerful enough to calculate the exact XC functional 
on systems of 100 atoms or more|25|. We will see how the 
exact functional deals with strongly correlated insulators 
and which new approximations are working for the right 
reasons. 

For more than a decade, researchers have been per- 
forming DFT calculations of molecular conductance, cal- 
culating the current in response to a bias applied to a 
molecule caught between two metal leads [50] . This prob- 
lem is prototypically difficult for present electronic struc- 
ture methods [5T]. This is a steady state situation, not a 
ground-state one. Model Hamiltonians typically used to 
study this kind of physics are insufficiently accurate, as 
several hundred atoms must be treated to achieve chem- 
ical realism [52]. A standard approach, combining the 
Landauer formula [53] with ground-state DFT and often 
called non-equilibrium Green's functions (NEGF) [M] , is 
typically used, and avoids adjustable parameters. But 
especially in the case of organic molecules, there is every 
reason to believe that, for the reasons given above about 
bulk gaps, our standard approximations are simply wildly 
inaccurate for this problem. Recent calculations [55] [55] 
including corrections based on non-DFT many-body cor- 
rections [57], and simple exact results for simple models 
[55] . are strong evidence that, at least for weak bias, the 
standard approach should yield accurate currents, once 
better DFT approximations are used [59] . 



IV. BACKGROUND: EXACT THEORY 

A KS calculation appears very similar to a Hartree- 
Fock calculation, but there's a crucial difference. The 
effective potential is defined to be one which makes the 
one-particle density p(r) be the exact density of the sys- 
tem. By virture of the HK theorem [23j [60], one can 
write 



E = T S + U + V nuc + E xc [p] 



(1) 



where T s is the energy of the KS orbitals, U is the Hartree 
(a.k.a. Coulomb) energy, V nuc is the attraction to the 
nuclei, and E xc [p] is everything else, defined to make the 
above exact. Then one finds the minimizing orbitals are 
given by the celebrated KS equations [24J : 



and 



,3 , MO 



(2) 



v s (r) = v nuc (r) + / d 6 r + v xc [p](r) (3) 

J \r- r\ 

and u xc (r) = SE XC /Sp(r). Thus, if the XC energy is 
known as a functional of the density, these form a closed 
set of self-consistent equations yielding the exact answer 
to the electronic structure problem, without ever calcu- 
lating the electron-electron repulsion directly! 
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FIG. 2. Exact radial density (red) and potentials (nuclear is 
blue, KS is pink) for the He atom in atomic units (Bohr radii, 
0.529A.). Two non-interacting fake electrons, doubly occu- 
pying the Is orbital of the pink potential, produce the exact 
ground-state density shown above. I thank Cyrus Umrigar 
for this density [61], 



In Fig. [2] I show the exact KS potential for a He atom. 
The blue line in the bottom half is v nuc (r) — —2/r in 
atomic units, the attraction of the electrons to the nu- 
cleus. In the top half, the red line indicates the exact ra- 
dial density of the He atom, found by using sophisticated 
wavefunction techniques to directly solve the Schrbdinger 
equation |62| . But the pink dashed line in the bottom is 
then the unique potential experienced by two fictitious 
non-interacting electrons that makes them have the ex- 
act ground-state density, i.e. 0i s (r) = \J p(r) jl. Every 
practical KS DFT calculation approximates v s (r), and 
never calculates the true many-electron wavefunction. Q 

Several relevant points about these equations include: 

★ Failures of DFT are due to failures of approximations. 
The exact E xc [p] yields E and p(r) exactly. 

■k Ground-state DFT calculations, even with the exact 
XC, yield only E and p(r) and any property that can be 
extracted from them (such as the ionization potential). 

★ The KS energies and orbitals replace those of HF for 
understanding chemical reactivity [63] , even though they 
appear as mere constructs for p(r). We also now know 
that orbital energy differences approximate optical tran- 
sitions, via TDDFT (Sec|W|. 

★ For most applications, we care only about £ as a func- 
tion of the nuclear coordinates. This determines bond 
lengths, angles, vibrational frequencies, all reaction en- 
ergies, transition-state barriers, etc. , i.e., most of the 
properties of interest about a molecule or a material. 



1 If you wish to annoy and confuse a traditional quantum chemist, 
ask "How much correlation is there in the KS wavefunction?" 
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V. CHALLENGE: LARGER SYSTEMS 

Although KS-DFT is just about the quickest way to get 
useful quantum calculations of electronic structure, its 
still too slow for a crucial application: MD simulations 
of liquids. There is tremendous interest in performing 
calculations that mix KS-DFT with other methods. For 
example, molecular mechanics with classical force fields 
can routinely handle millions of atoms, but making and 
breaking bonds is usually unreliable. On the other hand, 
the field of ab initio molecular dynamics (AIMD), using 
the Car-Parrinello method [64] of combining DFT with 
MD, is thriving, but limited to a few hundred atoms per 
simulation. Thus there is tremendous desire to perform 
mixed QM/MM simulations, in which the chemically ac- 
tive part of a larger system is treated in DFT. Many 
methods have been proposed, perhaps the most popular 
being ONIOM [65 , but many questions remain, espe- 
cially when the interface between the classical and DFT 
regions involve covalent bonds. This is a very active area 
of research. 

A related theme is that of orbital-free DFT, which 
looks backward to the original form of a pure DFT, as 
approximated by TF theory, and shown to be formally 
exact by the HK theorem. In a modern context, this 
means producing a sufficiently accurate approximation 
(much more accurate than TF theory) for the KS kinetic 
energy (T s in the language of Sec. IV), thereby avoid- 
ing the need to solve the KS equations. If this could be 
done, then all DFT calculations would run much more 
quickly, and AIMD could be performed on much larger 
systems. This has been an active area of research for 
many decades , with various empirical approaches be- 
ing tried for such approximations. Most (if not all) at- 
tempts use reasoning similar to that used successfully to 
approximate XC contributions, but there is nothing to 
suggest (and much evidence to the contrary) that such 
methods can produce a general purpose approximation 
sufficiently accurate for this purpose. Possibly a lower 
degree of accuracy could be acceptable for the larger sys- 
tem, and then a traditional KS treatment applied to the 
chemically active region. This embedding idea has been 
tried, especially when the system is weakly bonded to the 
environment |66| 167] . 

Many of these ideas are unified in a new approach to 
the issue of subsystems within electronic structure, called 
partition DFT (PDFT). Originally developed by Cohen 
and Wasserman to resolve difficult issues in chemical re- 
activity theory 68 , this formalism unifies several distinct 
concepts in chemistry, including atoms in molecules, ef- 
fective charges, chemical identity, localization of bonds, 
etc. To date, only calculations on model systems [ST)] and 
diatomic molecules [70] have appeared, but interest is 
growing rapidly [71] An interesting challenge is to adapt 
PDFT to include external electric fields, as it seems ide- 
ally suited to the molecular electronics problems of Sec. 

inn 



no simple rule for reliability 
no systematic route to improvement 
decades between each generation 
full of arcane insider jargon 
too many approximations to choose from 
can only be learned from a DFT guru 



TABLE I. List of things users despise about DFT calculations. 
Please rank in order of induced frustration. The extreme left 
column indicates the users response to these 'features', the 
right denotes that of developers. 



VI. BACKGROUND: THE USERS' LAMENT 

But the popular success of DFT has bred its own set 
of problems. As codes become faster and easier to use, 
DFT is applied to a huge number of situations. Fig. [T] 
includes applications to protein folding [72] , astrophysics 
[73] . dyes [74], and dirt [75], to name just a few. As ex- 
perience builds with each given functional, the accuracy 
and reliability comes to be known, as well as the quali- 
tative failures. For example, it was very early recognized 
that standard approximations do not yield long-range 
dispersion [76] . The local nature of the standard approx- 
imations implies an exponential decay of the interaction. 
Progress in this area only came recently, and is described 
in Sec. [X] Another example is that of anions, which are 
technically unbound with standard approximations [77] . 
and yet for which accurate results can still be obtained 
[78] . Only recently has this puzzle been addressed and a 
solution proposed [75] , 

Naturally, users would like an all-purpose tool that 
provides answers of a prescribed quality in all situations. 
Present DFT calculations are a far cry from this. Table 
I shows a list of features that most users will unfortu- 
nately recognize. Throughout its history, DFT has pro- 
vided approximations that work for some problems and 
fail for others, in largely mysterious ways. At any given 
moment, the most popular approximations fail for the 
most interesting systems, such as the moderately corre- 
lated oxides of Sec. IIII1 



VII. PROGRESS: EXCITATIONS 

Because an excited-state density does not uniquely de- 
termine the potential [80], there is no general analog of 
HK for excited states. Many different ways to approach 
the calculation of excited states in DFT have developed 
over the years, including ensembles [EH E2] ASCF [55] . 
min-max principles 84] , and others. However, since the 
mid-90's, TDDFT has become extremely popular [85] , 
Because the methodology is technically very similar to 
that of TDHF, TDDFT was very rapidly implemented 
in quantum chemical codes such as Turbomole [86], and 
is now a standard part of any code. A useful reference 
and pedagogical tool is [87 . 
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TDDFT is based on a similar (but distinct) theorem 
to HK. he Runge-Gross theorem [88] establishes that, in 
a time-dependent quantum problem, all observables are 
functionals of the time-dependent density (under certain 
conditions). Then, considering the linear-response of a 
molecule to a time-dependent electric field, one finds sim- 
ple formulas that correct the KS eigenvalue differences 
into the (in principle, exact) optical excitations of the 
system |89| . but requiring (of course) another unknown 
functional, called the XC kernel [90]. All practical calcu- 
lations employ the adiabatic approximation, and almost 
all use the same approximation for both the ground-state 
and the TDDFT calculation [5T]. 

Practical TDDFT often produces good excitation spec- 
tra. Typical errors in individual energies are higher than 
for the ground-state (0.4 eV), but properties are excellent 
(bond lengths, dipole moments, etc.), and even roughly 
accurate spectra can be sufficient to identify the domi- 
nant excitations in the optical spectra of large molecules 
[9"T] . About f0% of all DFT calculations now include 
TDDFT as well, as its computational cost is not many 
times more than a single calculation of a ground-state 
energy. The procedure can even be extended to contin- 
uum states [92J to accurately predict electron-atom elas- 
tic scattering. 

But, just as with ground-state DFT, along with suc- 
cess come challenges. The standard functionals are in- 
accurate for charge transfer excitations, leaving them so 
low so that they contaminate other parts of the spec- 
trum. Double excitations are also excluded by the adia- 
batic approximation. Calculations of the optical response 
of solids look very similar to RPA (see below) if using 
standard approximations, which are insufficiently non- 
local in the case of a solid. For applications beyond lin- 
ear response (strong fields), the standard approximations 
cannot be used because of their poor quality potentials. 
Furthermore, sometimes the quantity of interest is not 
a simple functional of the time-dependent one-electron 
density, such as the double ionization probability 93 . 

But all these are open areas of research and progress is 
being made continuously. Many papers have been pub- 
lished on the charge transfer issue, and several range- 
separated schemes seem to handle these excitations well 
[M] , The basic form of the kernel needed to handle dou- 
ble excitations was deduced years ago [55] , and has been 
generalized and systematized successfully [95] . With 
much effort, the nature of the error for solids was de- 
duced [97] . and a new approximation to the kernel (the 
nanoquanta kernel) can be found from many-body per- 
turbation theory [55] . In a very recent development, a 
bootstrap approximation for the kernel in terms of the 
dielectric function appears to produce accurate excitons 
in all but the largest-gap insulators [99] , 



VIII. BACKGROUND: FIRST PRINCIPLES OR 
UNPRINCIPLED? 

Over the decades since the introduction of the 
Schrodinger equation, many excellent methods have been 
developed for directly solving the electronic structure 
problem, including configuration interaction, Moller- 
Plessset perturbation theory, the coupled-cluster expan- 
sion, and quantum Monte Carlo |I00j . Such methods 
rarely suffer from any of the difficulties listed in Table I. 
But because of the coupling between coordinates in the 
many-electron Schrodinger equation, the computational 
cost of such methods is usually significantly higher than 
that of DFT. Loosely, the more accurate the method, 
the more rapidly the cost rises with number of atoms. 
Thus, without some algorithmic breakthrough, DFT will 
always allow more atoms (often by a factor of f 0) to be 
treated, no matter how fast our processors get, or how 
many we have. 

This does not make direct solutions obsolete. They 
provide crucial benchmarks for testing approximate func- 
tionals for smaller systems, and give crucial insight into 
the nature of errors, both quantitative and qualitative. 
In chemistry, it is traditional to refer to standard ap- 
proaches as ab inito, while DFT is regarded as empirical. 
Because solid-state calculations are more demanding, for 
many decades DFT was the only possible approach. Thus 
DFT calculations are referred to as ab initio in solid-state 
physics and materials science. This is why there is a 
solid-state code ABINIT [TUT] which performs only DFT 
calculations. 

A sore point is whether or not approximate DFT 
should be called empirical. Even if an approximate func- 
tional includes parameters that have been fit to some 
data set, once the final form has been written down, that 
approximation can be applied to every possible electronic 
structure problem, without adjusting parameters for each 
specific calculation. Thus DFT, with a fixed approximate 
functional, is still first principles, in the sense that the 
user only chooses the atoms, and the computer predicts 
all properties of the molecule or solid. 

As mentioned above, the first approximation was LDA 
|24| , and the formula for this is determined by properties 
of the uniform electron gas. No-one disputes that DFT 
with LDA is non-empirical [102, I03J. But even just the 
next step up Jacob's ladder [35] of functional sophisti- 
cation, the GGA has no unique form. There are two 
major schools of thought here. Purists like to use exact 
conditions of quantum mechanics to derive the param- 
eters in their approximate functionals, and so claim to 
be non-empirical. This school has been championed by 
John Perdew, with a lifetime of very successful approx- 
imations [5] 11031 1104] . On the other hand, pragmatists 



2 When teaching, I explain that DFT is some acronym for unreli- 
able, while ab initio is Latin for too expensive. 
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like Axel Becke and Bob Parr have allowed one or two 
parameters to be fit to specific systems, such as in B88 for 
exchange [7] and LYP for correlation [3]. Such approx- 
imations have been based on sound physical reasoning 
underlying the structure of the approximation. I have 
even had the pleasure of deriving some of these parame- 
ters much later |105l 1106] . By fitting, one usually finds 
higher accuracy for systems similar to those fitted (often 
by a factor of 2), but greater inaccuracies far away. For 
example, LYP correlation [8] works very well as part of 
B3LYP in chemistry, but fails badly for bulk metals. The 
PBE approximation [9 works passably well for many ma- 
terials purposes, but can be a factor of 2 or more worse 
than BLYP for dissociation energies. 

A simple example of a first-principles approach is given 
by the B3LYP approximation (and its materials counter- 
part, PBEO [107] ). The crucial part of a hybrid, which 
mixes HF exchange with a GGA, is the fraction of exact 
exchange, a, which is about 20%. This was fixed once 
and for all in the definition of the functional, and all the 
thousands of papers using B3LYP in Fig. [l]are gather- 
ing information on the same functional. The amount of 
mixing can be rationally related to other calculations of 
atomization energies [108] . If authors adjust the amount 
of mixing to improve their results for some system or 
property, this is not first principles. 

Many of these tensions have been highlighted in recent 
years by the Minnesota functionals developed by Truhlar 
and co-workers |109[ 1110] . These use the same basic in- 
gredients as the standard approximations, but optimize 
performance on a training set of energies by fitting up to 
several dozen parameters. They often produce more ac- 
curate results on systems close to those trained on, but 
yield little insight into functional construction and de- 
velopment. I consider these first-principles, because the 
functionals contain no parameters adjusted to the system 
being calculated. But I also avoid altogether the e-word, 
because such functionals illustrate the uselessness of that 
word in the context of DFT approximations. [^] 

Fig. [3] shows a random selection of some of the hun- 
dreds of XC approximations that have been suggested 
over the years. Many of these are in popular codes, some 
of which even allow you to design your own functional. 
Clearly, when calculating a property to compare with ex- 
periment, one could keep trying functionals until agree- 
ment with the measured value is reached. Not only is this 
contrary to the entire spirit of DFT, it is certainly not 
first-principles, and is the worst form of empiricism. The 
literature today is rife with such calculations, and the 
existence of so many approximations, with so little guid- 
ance, is bringing the field into disrepute. Users should 
stick to the standard functionals (as most do, according 
to Fig. [lj, or explain very carefully why not. 
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3 Interestingly, I have recently co-authored an approximation with 
about 10 'empirical' parameters [111] , Don Truhlar, eat your 
heart out! 



FIG. 3. The alphabet soup of approximate functionals avail- 
able in a code near you. Thanks to Peter Elliott for the 
artwork. 



IX. CHALLENGE: DIGGING DEEPER 

As we have seen, the practice of modern DFT suffers 
a lack of detailed understanding of how to approximate 
functionals. We begin from local approximations, and 
then create more accurate, sophisticated versions. Based 
on insight and intuition, combining either exact condi- 
tions from quantum mechanics or fitting specific sys- 
tems, we make progress, but only very slowly, and rarely 
without ambiguity. A formally exact theory exists (Sec. 
IV), but provides only limited guidance about approxi- 
mations. 

I believe that a fundamental principle underlies the 
success of DFT, which is that local approximations are 
a peculiar type of semiclassical approximation to the 
many-electron problem [112] . For the last 6 years, with 
both my group and many collaborators, I have been try- 
ing to uncover this connection, and make use of it. The 
underlying math is very challenging, and some must be 
invented. 

Basic quantum textbooks have separate chapters on 
perturbation theory and semiclassical approximations, 
but never relate the two [113] . All modern many-body 
methods have their roots in treating the interaction as a 
perturbation, since we find solving non-interacting prob- 
lems relatively easy. But such treatments ignore the 
fact, proven by Lieb more than a quarter of a century 
ago |114j . that TF theory (see Sec. [TTJ) becomes rela- 
tively exact for neutral atoms as Z — > 00. As detailed 
by Schwinger [115] and others, this is the semiclassical 
limit [1161 1117] mentioned above. Dramatic confirmation 
of this fact is that TF also appears to yield the exact ion- 
ization potential of atoms (averaged across a row) in this 
limit, and that LDA-X recovers the HF result, including 
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the oscillations across a row |118j . Most of our results so 
far have been confined to Id systems |116) . but this in- 
sight lead to the restoration of the gradient expansion in 
PBEsol, which cures many of the PBE problems with lat- 
tice parameters |119j . and a derivation of the parameter 
in the B88 functional [TUB] . 



X. PROGRESS: WEAK INTERACTIONS 

The ability to treat van der Waals is a recent (and 
ongoing) success story for DFT. In the 1990's, it was 
well-known that the standard functionals could not yield 
correct long-range dispersion forces [76], i.e., their bind- 
ing energy curves decay exponentially (with density over- 
lap) instead of —Cq/R 6 , where R is the separation and 
Cq is the van der Waals coefficient, determined from the 
the frequency-dependent polarizabilities of the fragments 
[120) . Because this excluded such a huge number of im- 
portant systems and properties (such as DNA, physisorp- 
tion on surfaces, most biochemistry, etc.), there were al- 
ways adhoc methods for adding back in dispersion, using 
pair potentials between atoms [121] . 

Over about 20 years of research and many papers, 
the late David Langreth with Bengt Lundqvist [122] 
and many other collaborators, developed an approximate 
non-local ground-state density functional, call it LL, that 
has the right decay behavior and reasonably accurately 
captures these effects |123| . LL is entirely non-empirical, 
using results from the uniform gas and interactions be- 
tween slabs of that gas to find such a form from first 
principles. While the initial implementation of this func- 
tional was computationally expensive, a recent algorithm 
of Soler |124j made it much faster, so much so that its 
cost is negligible beyond about 100 atoms. This led to im- 
mediate implementation in many codes worldwide, and 
there is now a plethora of calculations with LL [125] . 

Simultaneous with this development, in quantum 
chemistry, Grimme |126j developed his DFT-D method- 
ology that provides an empirical correction to DFT re- 
sults in a systematic and accurate fashion. The results 
for small molecules in the S22 data set |127j are ex- 
tremely good. DFT-D is much more accurate for these 
systems than LL, but LL can be applied to all matter (ex- 
cept -possibly- metals), including situations where pair- 
potentials cannot work. A less empirical alternative to 
Grimme has been proposed by Tkatchenko and Scheffler 
[128] , which produces a scheme for calculating an additive 
correction for any functional and has only slight empiri- 
cism. This has recently been extended to include even 
metals [129] . Over the next five years, one (or possibly 
two) of these schemes is likely to become the standard 
method for including weak interactions in DFT. 



XI. CHALLENGE: TIME FOR A PARADIGM 
SHIFT? 

The KS equations, combined with the local density 
approximation, were a reformulation of the electronic 
structure problem relative to TF theory. By producing a 
more demanding computational algorithm while lessen- 
ing the fraction of the total energy that needs approxi- 
mating, a great leap forward in accuracy and reliability 
was achieved. But that was back in 1965. Perhaps we 
are at the end of that road in terms of useful approxima- 
tions, and what is needed now is a new paradigm which 
begins from different starting point. 

Optimized effective potential (OEP): [1301 1131] At 
some point, exact exchange (loosely, evaluating an 
orbital-dependent functional in the KS scheme) seemed 
like a strong candidate, because it allowed exchange to be 
evaluated exactly, instead of being approximated. This 
cures a multitude of problems with local and semilocal 
approximations, such as better orbital energies, avoid- 
ing self-interaction error, and better approximations to 
the derivative discontinuity. However, the technology for 
solving the OEP equations efficiently has existed for at 
least a decade [132H134] , but no general-purpose correla- 
tion approximation has been found that works well with 
exact exchange in all situations. Even worse, there are 
questions about how well-defined the equations are in a 
basis [135] . 

Random phase approximation (RPA): A more recent 
development |136j . driven primarily by improvements in 
hardware and algorithms, is the ability to solve the RPA 
equations efficiently for systems up to about a hundred 
atoms. These can now be done faster than a conventional 
Hartree-Fock calculation. Moreover, bare RPA has many 
excellent features, including exact exchange and a van der 
Waals contribution, but also has problems with atomiza- 
tion energies. RPA is now implemented in both materials 
[137] and chemical codes [25]. RPA makes a promising 
candidate for a new baseline calculation, to which further 
inexpensive approximations can be added, because it in- 
corporates more traditional theory (here, coupled-cluster 
|138j ) at low computational cost. 

Density matrix functional theory (DMFT): A darling 
of the chemistry community, and of increasing interest 
in physics, is density-matrix functional theory (DMFT) 
[139] . in which impressive results have been gotten over 
the years. The formalism is well-founded, using simi- 
lar variational principles as DFT, due to Gilbert |140j . 
and a sequence of successive approximations and re- 
finements have produced excellent results at equilibrium 
bond lengths and for total energies of closed-shell sys- 
tems [141] . Even gaps of insulators seem to come out well 
[142] . But the methods are not yet developed for general 
purpose, as open shell systems and size-consistency re- 
main issues. If this methodology ever does become pop- 
ular, it would represent a true paradigm shift, as it does 
not even use KS equations. But, for this reason, it is dif- 
ficult to see how many of the impressive results of DFT 
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approximations could be retained. 



XII. THE FUTURE 

So, where does this leave us? It is clearly both the 
best and worst of times for DFT. More calculations, both 
good and bad, are being performed than ever. One of the 
most frequently asked questions of developers of tradi- 
tional approaches to electronic structure is: "When will 
DFT go away?". Judging from Fig. [T] no time soon. 



Although based on exact theorems, as shown in Fig. [2] 
these theorems give no simple prescription for construct- 
ing approximations. This leads to the many frustrations 
of the now manifold users listed in Table [U Without such 
guidance, the swarm of available approximations of Fig. 
[3] will continue to evolve and reproduce, ultimately un- 
dermining the entire field. Let us hope that some of the 
many excellent ideas being developed by the community 
will come to fruition before that happens. 

I am grateful for the support of NSF under grant no. 
CHE-1112442, and the assistance of past and present stu- 
dents with the preparation of this perspective. 
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